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ABSTRACT 

Methods of summation of power series relevant to applications in quantum theory are 
reviewed, with particular attention to expansions in powers of the coupling constant and 
in inverse powers of an energy variable. Alternatives to the Borel summation method 
are considered and their relevance to different physical situations is discussed. Empha- 
sis is placed on quantum chromodynamics. Applications of the renormalon language to 
perturbation expansions (resummation of bubble chains) in various QCD processes are 
reported and the importance of observing the full renormalization-group invariance in 
predicting observables is emphasized. News in applications of the Borel-plane formalism 
to phenomenology are conveyed. The properties of the operator-product expansion along 
different rays in the complex plane are examined and the problem is studied how the 
remainder after subtraction of the first n terms depends on the distance from euclidean 
region. Estimates of the remainder are obtained and their strong dependence on the na- 
ture of the discontinuity along the cut is shown. Relevance of this subject to calculations 
of various QCD effects is discussed. 
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1.1 Perturbation series 

A typical situation in physics is want of exact solutions. To find a good approximation of 
a problem, we have mostly first to simplify it by neglecting a number of effects, thereby 
facilitating its solution but, simultaneously, endangering its physical relevance. It is a 
general experience that exact equations can mostly be solved only approximately. 

Perturbation theory (PT) is based on the idea of expressing the solution f(g) to a 
problem (equation) in the form of the power series 

oo 

E fn9 n , (1) 
n=0 

where g, the perturbative parameter, is considered to be a small numerical quantity 
describing, possibly, a physically important effect. 

The first question to ask is how f(g) is related to the series ([!]). To approach it, we 
have to see (i) what is the precise meaning of (p]) and, (ii), in what sense the series ([[]) 
is assigned to /(<?), the searched solution to the problem. For this to discuss, we have to 
consider complex values of g. 

Simplest is the situation when f(g) is holomorphic inside a circle of radius p > centred 
at the origin; then (^) is the Taylor expansion of f(g) and determines f(g) uniquely, f(g) 
being equal to ([TJ) for all \g\ < p. The /„ are given by the derivatives of f(g) at the origin. 

If, however, p = 0, the equality 

oo 

/(<?) = E fn9 n , (2) 
n=0 

being valid only for g = 0, is not relevant for applications. On the other hand, if we 
assume, instead of (Q), an asymptotic relation between f(g) and ([!]), 

oo 

/(<?) ~ E fn9 H (3) 
n=0 

for g — > 0, there is always an f(g) satisfying (|3]) for any set {/„} |IJ. Further, in contrast 
to (0), the relation (^) does not determine the function f(g) uniquely, even when all the 
f n are explicitly given and the set of rays approaching the point g = is specified. Even 
then a whole class of functions f(g) is defined by ([3]). The situation changes if additional 
restrictions are imposed on the remainder f(g) — J2n=o fn9 n ] then there may be only one 
function satisfying the requirements, and we can call it the sum of the series. 
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Such additional restrictions as well as the assumptions (|2|), ([|) themselves should be 
motivated physically. Equation (0) is related to the question of convergence: 

• /l/ (Convergence): Is the series convergent, in the sense that the limit of 

JV 

/(*)(<?) = £/n<f (4) 
n=0 

for iV — > oo, with g fixed, exists? For what values of g is this the case? 

The asymptoticity property (|3|), on the other hand, is a statement on smoothness: 

• /2/ (Smoothness): Is the expression 

N 

\f(g)-Y.fng n \ (5) 

n=0 

bounded, in what domain of (N, g) (with g complex and N a nonnegative integer) 
and what is the bounding function? 

Finally, there is the question of uniqueness: 

• /3/ (Uniqueness): Can the sequence {f n }, n = 1,2,..., uniquely determine the 
expanded function f(g) and under what conditions? 

Perturbation theory yields, at least in principle, explicitly all the coefficients /„ of the 
expansion ([!]). The knowledge of the large-order behaviour of the f n ,n = 1,2,..., tells 
us whether the series (|1|) is convergent or not, but what we really need is to answer the 
question /3/, i.e., that of a unique determination of f(g). If f(g) is holomorphic at the 
origin, ([]]) is convergent and f(g) is uniquely determined from the coefficients /„ by means 
of (0). There was a universal belief till the early 1950's that a sound physical theory should 
not use a divergent perturbation expansion or a function singular at the origin.Q 

In 1952, F.J. Dyson || argued that in QED there is a singularity at the origin of the 

complex coupling-constant plane. (The conjecture that such a singularity makes the series 

divergent was questioned later, see PJ and references therein.) 

1 This opinion survives among many physicists till the present time. There was a similar mistrust 
of the use of divergent series also among mathematicians, but much earlier. J. Basdevant gives two 
interesting passages, by J. d'Alembert and by N.H. Abel. Other examples illustrating the evolution of 
opinions on divergent series are in ref. [GJ. 
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A wide class of divergent series can be given precise meaning by means of Borel sum- 
mation, provided that certain additional conditions are imposed on f(g). This method 
has found fruitful applications in quantum field theory, where its use goes back to the 
1950's. In 1951, for example, J. Schwinger || obtained a compact expression for the term 
in the Lagrangian corresponding to the QED vacuum polarization by an external constant 
magnetic field H. Then V. Ogievetsky used the (divergent) perturbative expansion of 
this term to show that its Borel sum coincides with Schwinger's result. Indeed, the coeffi- 
cients of its perturbative expansion in powers of e 2 behave as (2n — 3)!(— H 2 / (7r 2 m A )) n at 
large n, where e and m is the charge and the mass of the electron respectively. Ogievetsky 
found the Borel sum in the form 

roc „ 

/ e- mt a(t,e)dt (6) 
Jo 

with 

a(t,e) = --JL_[efcff cot{etH) - 1 - (etH) 2 /3), (7) 

which coincides with the compact expression obtained by Schwinger ||. This important 
result shows that a divergent perturbative expansion does not signal an inconsistency in 
a theory; it also shows that there are special - but realistic - cases of Borel summability 
in QED, although general considerations indicate Borel non-summability (see [|7|, ||, [[(J, 
and Table 1 and a discussion in section 2 of the present paper). 

Gradually, the Borel summation techniques became widely adopted in quantum theory. 
It even became common to judge theories according to whether their perturbation expan- 
sions are Borel summable or not; in this sense, the formal criterion of Borel summability 
became successor to the criterion of convergence. Being rather technical assumptions, 
however, neither convergence nor Borel summability can serve as criteria of consistency of 
a theory. One might find it more natural to set another criterion of consistency: that the 
transition from the perturbed to the unperturbed state should be continuous and smooth 
in the coupling parameter. But even this is not necessary: Nature is not obliged to be 
smooth in our parameters. It is a special step to assume that the transition from the 
world with interaction, g > 0, to that without it (g = 0) is smooth. This is why the 
relation (|3|) has nontrivial physical implications. 

Disregarding special field theories with interactions being interpreted as discontinuous 
perturbations in the sense that the theory does not return to the unperturbed theory 
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as the interaction coupling vanishes, the fact is that asymptotic expansions in powers 
of g have been universally adopted in quantum theory; perturbation theory does imply 
the assumption of a smooth transition from the state with interaction to that without 
it. What we then really need is to answer the question /3/ of uniqueness, which is 
independent both of convergence and of Borel summability. If the answer to /3/ is yes, 
i.e., if the /„ determine f(g) uniquely, the theory could be considered as selfsufficient, no 
matter whether the conditions of Borel summability are satisfied or not. If the answer is 
no, the theory has to be supplemented with some additional input. 

In the majority of situations in physics, f(g) has a singularity at the origin; note that 
smoothness at a point does not exclude a singularity at this point. Then the perturbation 
expansion is not a Taylor series and, if considered as an asymptotic expansion, it needs a 
supplement for uniqueness to be reached. Needless to say, this supplement should be of 
physical nature, while no formal argument (e.g., that the method yields a finite value of 
a divergent series) can fill the gap where physics is missing. 

1.2 Large-order behaviour 

The problem of the high-order behaviour of the perturbation-expansion coefficients in field 
theory calculations has aroused new interest in the last decade. Particular attention has 
been paid to power corrections to QCD predictions for hard scattering processes. One 
can point out two reasons for this growing interest: 

• The theoretical problem of how physical observables can be reconstructed from their 
(often divergent) power expansions. 

• The pragmatic need to assess the usefulness of performing the extensive evaluations 
of multi-loop Feynman diagrams in QCD. Much effort has been devoted to the 
computation of higher-order QCD perturbative corrections; in some cases third- 
order approximations are known, and we now seem to be at the border of what can 
be carried out analytically or numerically in high-order perturbative calculations. 
So, why next order? If the series is divergent, the next order may represent no 
improvement with respect to the lower-order result. On the contrary, at a certain 
order it will lead to deterioration. 
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The question of the possible divergence of asymptotic perturbation expansions, which 
goes back to the original argument of Dyson [Q, still defies definite and clear answer. 
Notice in this connection that a series that is asymptotic to a function singular at the 
origin need not be divergent. Indeed, for any series ([!]) convergent within a circle \g\ < p 
there are functions f(g) which are singular at g — and have (|1|) as asymptotic series. The 
traditional wisdom that the divergence of perturbation expansions is directly related to the 
discontinuities of Green's functions is not generally true. If the series is asymptotic, the 
knowledge of the large-order behaviour of perturbation-theory coefficients is insufficient for 
the determination of the function expanded and one needs much more detailed information 
in the vicinity of the point of expansion in order to establish the relation between the 
expansion and the function expanded. 

In fact most of the information on the analytic properties of relevant Green's functions 
near g — || [Tl| may actually be invisible in perturbation expansions themselves. More- 
over, the summability of a perturbation expansion is not the same as the possibility of 
recovering from such series the corresponding full physical quantity. As the summability 
of QCD perturbation expansions is determined by the large-order behaviour of their co- 
efficients, we primarily need to know which kind of information is necessary to explicitly 
obtain the sum of the series. 

If a theory provides, in addition to the values of the expansion coefficients /„, enough 
supplementary information, a power series can be summed to a finite number, which is 
not arbitrary but reflects the additional piece of information used. One has to give it such 
a form that it could be used to reduce (or even remove) the ambiguity connected with the 
asymptoticity of the expansion. The task of a summation method then is not to assign 
arbitrarily an ad hoc value (whose only merit might just be that it is finite) to your infinite 
or ambiguous series, nor to fill the gap in physics (which your theory possibly possesses) 
with some mathematical cement, but to offer you a service: to translate your physical 
idea (or that of the theory you work with) into a mathematical condition to reduce the 
ambiguity of the formal power series. 

As a typical example of this procedure the role of general principles in the S-matrix 
theory can be mentioned. The S"-matrix techniques were built from mathematically for- 
mulated, but physical principles: analyticity as a mathematical formulation of microscopic 
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causality, unitarity as probability conservation, crossing symmetry as a general property 
of Feynman diagrams, and the position and character of the cut in the complex energy 
plane as a description of the properties of the particle spectrum. A simultaneous use of 
the analyticity properties and of the concepts of field theory was carried out in the QCD 
Sum Rules program ||12| , which met with tremendous phenomenological success. Note 
that the foundation of this approach rests on the non-trivial fact [O, O] that the princi- 



ples of S-matrix theory are rigorously true in QCD: it follows that for gauge theories, 
if confinement is imposed, the analytic properties of physical amplitudes are the same as 
those obtained from an effective theory involving only the composite (physical) fields. 

It should be considered as very fortunate that, simultaneously, analyticity plays a 
crucial role also as a mathematical condition reducing the ambiguity of asymptotic series. 
In Section 2 of the present paper, we discuss in detail the interplay between large-order 
behaviour of a series (as listed in Table 1) and the analyticity properties of the function 
expanded; it turns out that a balance between these two concepts is needed for a unique 
determination of f(z) from @, in the sense that if more analyticity of f(z) is available, 
one can afford a more violent behaviour of the a n , and vice versa. 

In Section 3 we focus on some practical aspects of the operator-product expansion, 
in particular on the problem of how the remainder after subtraction of the first n terms 
from the function expanded depends on the distance from euclidean region, provided that 
an estimate on the remainder in euclidean region is known. Sections 4 and 5 are de- 
voted to singularities in the Borel plane and applications of the Borel plane formalism to 
phenomenology. We discuss the progress in understanding the role of the singularities in 
the Borel plane for nonperturbative effects, in studying hadrons containing heavy quarks, 
then the different methods of calculation of the coupling constant a s from existing ex- 
perimental data, the role of operator product expansion in practical calculations and the 
progress in the resummation of bubble chains. Section 6 is devoted to problems related 
to renormalization ambiguities and their impact on phenomenology. Section 7 contains 
concluding remarks. Some parts of the present paper represent an extended and amended 
version of ref. |15 |. 

Below we give a survey of the large-order behaviour of perturbative series in different 
theories. Bender and Wu pi examined perturbation theory for one-dimensional quantum 
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mechanics with a polynomial potential and they obtained the large-order behaviour of the 
perturbation coefficients. Lipatov |Tj| obtained the same results for massless renormaliz- 
able scalar field theories. 

Brezin, Le Guillou and Zinn- Justin ||19|| , applying independently the same method to 
anharmonic oscillations in quantum mechanics, were able to rederive and generalize the 
results obtained by Bender and Wu. 

In QED, after the pioneering work of Dyson a number of papers appeared discussing 
the analyticity properties of the Green functions (some of them are cited in refs. H and 



||20|). The growth like n\ (where n is the order of approximation) has two sources: 

1. The number of diagrams grows like n\, each diagram giving a contribution of the 
order of 1. 

2. There are types of diagrams for which the amplitude itself grows like n\ f2l| . 

A survey of the large-order behaviour of expansion coefficients in some typical theories 
and models is given in the Table 1. As subtle cancellations among higher-order graphs 
may occur, the expressions in the third column of Table 1 may sometimes give an upper 
bound rather than the actual high-order behaviour of the coefficients. 

Table 1 is intended for first information and should not be used for systematic anal- 
yses because some important conditions or restrictions could not be mentioned. A brief 
explanation of its use is given below. Notation is not consistent, being mostly taken from 
the papers quoted; the reader is referred to them for explanation. References in the table 
as well as in the text are often made not to the original papers but rather to more recent 
papers or reviews, in which the subject is explained in modern context. As a result, a 
number of valuable papers are not quoted. I apologize for this. For details I refer to the 
papers cited and references therein, and to the anthology by Le Guillou and Zinn- Justin 



20| , which contains a list of references up to 1990. 

A few words on the explanation of the symbols used in the table. In the item 1, 
f(z) is a real function analytic in the complex z plane cut along [— oo,0] and equal to 
J oo dz'lmf(z')/(z' — z). In item 2, the E m (g) are, for g > 0, the eigenvalues of the 
Hamiltonian H = — |(d/dg) 2 + \q 2 + \gq 4 . 

Eckmann, Magnen and Seneor pSI consider the two-dimensional Euclidean boson field 
theories and give bounds on truncated Schwinger functions. They also prove that the 



S 



Theory 



Notation and references High-order behaviour 



1 Disp. relation 

2 Anharmonic 
oscillator 

3 Anharmonic 
oscillator, 4 

4 (p 4 in 2 and 

3 dimensions 

5 Instantons, 

D = 2m/(m-2) 

6 Field theories 
without fermions 

7 Field theories 
with fermions 

8 Yukawa theories 

d = 2 

9 QED 



10 QCD in 2 dim. 

11 QCD 

12 Bosonic strings 

Table 1: High-order 

(see also explanations 



Im/(z) ~ z- b e- a ' z 
E m (g) = J2T=o E m ,n9 n 



Hg) ~ f-oo ei ~ x /2 ~ 9X /4) dx 



25 



26 1 



m > 6, |JTJ 
m = 4, |] 

^ = Er=o2f ) (-eV(^)) fc 



0,111- 



li 

(1st IR renormalon) 

01 



46 1 



see [44], [45| for discussion 
is # of handles, p2fl 



/„~(-fl) n r(n + 6) 
£ m ,n ~ (-3/4)"n! 

I n ~ (-l) n 4""(n - 1)! 

Analyticity for complex g, 
Borel summability proved 

C n (m) ~ (_^)« e "(l-™/2) n (m+D)/2 

C n (4) ~ 2.75n 4 (n/(167r 2 e)) n 
Z< fc ) ~ A:!A; 6 "A- fc c„[l + 0(1/Jfe)] 

A n ~ r(n(d - 2)/d)K n n- a cos(2nn/d)A 

f n ~ (-l) n 5 rt / 2 r(n/2)D 

Z n ~ n- a A-"cos(2vm/rf)r(n(rf- 2)/d) 

d n ~ A~ n (ln n) n 

/ n ~(-l)"r((n + z/)/2)a"+^ 

<2 n ~ const(-/3 /2)«r(n + 1 - 

r n ~ const /3£n! 7l 2+/?i/^+72 

(^ 2 AT c /2) 2fc (-l) /£ - 1 (2fc)! 

/„ ~ a"n 7 n! 

rf n ~ const(-/3 /2)™r(n + 1 - 2/V$) 



behaviour of perturbation expansion coefficients 

in the text) 
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domain of analyticity and of the bounds obtained can be extended to an angle > it bisected 
by the positive real semiaxis, and prove the Borel summability of the power series of the 
Schwinger functions. In the case of the Euclidean 4 theory in 3 dimensions, Magnen 
and Seneor [2£| prove the stability of the free energy for complex values of the coupling 



constant and obtain, in analogy with the two-dimensional case, Borel summability of the 
perturbation series. 

In the instanton row (item 5), the C n (m) are the coefficients of the expansion of 
the Gell- Mann-Low function, with n being the perturbative order and m being related 
to the dimension D by D = 2m/ {m — 2). Further, both C n (m) and C n (4) contain an 
n-independent (but m-dependent) factor. 

The item 6 is related to scalar electrodynamics in four dimensions, with the action 



\Fl v + \\D^+\m^ + \\o\\iV. 



(8) 



The Z n are one-particle irreducible Euclidean Green's functions, with c n , b n and A being 
constants and A being equal to 12 approximately. 

For field theories with fermions (item 7), A(g) = J2 n A n g 2n is the Green function of the 
scalar field a with the coupling term ga^>^ for a Yukawa interaction in d > 2 dimensions 
H; R, a and A are computable constants. In ref. ]3~T|], the f n denote the main contribution 



to the QED perturbation theory asymptotics, where 5* is equal to 3 3 / 2 47r 3 . 

The d n in the item 9 describe (see the asymptotic behaviour of perturbation 



series due to the first infrared renormalon in QED. (In QCD, item 11, the result agrees 
with that obtained by relating infrared renormalons to nonperturbative corrections.) The 
d n are the coefficients in the perturbative expansion of the function Q 2 d/(dQ 2 ) P(Q 2 ) in 
powers of the coupling parameter a(Q), 

J2 d n a(Q) n+1 . (9) 

n=-l 

In QED, aP is equal to IT, the photon vacuum polarization, and gLi = — j3o, fio is the 
first coefficient in the expansion of the beta function, which is defined as /i 2 d/d/i 2 a(/i) = 
Pin) = E n =o PnQ-{[j) n+2 ■ In QCD, P = II is the electromagnetic current-current cor- 
relation function. The large-order behaviour in QCD is d n ~ const (—/?o/2)™r(n + 1 — 
2/3i//3o)(l+0(n)). Beneke and Zakharov [33| show that in QED (with massless fermions), 
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the coefficients d n have the same IR behaviour with j3 n taking their QED values. The UV 
renormalons are also present and produce the fixed-sign divergence of the perturbative 
series. They have been analyzed by Beneke and Smirnov who study classes of dia- 
grams with an arbitrary number of chains, and derive explicit formulae for leading and 
subleading divergence. To organize the diagrams in classes, the expansion parameter 1/Nf 
is used, where Nf is the number of fermion species; as a consequence, diagrams suppressed 
in the 1/Nf expansion are not suppressed for large n and, consequently, no finite order in 
the 1/Nf expansion provides the correct behaviour in n in the full theory. Table 1 shows 
the large-order behaviour of the vacuum polarization, r n being the coefficient of ia n+1 in 
the perturbative expansion and 72 = 99/(8iV|) . The authors discuss extension of the 
formalism to non-abelian gauge theories and expect a similar result. 

For QCD in 2 dimensions (with N c tending to infinity), Zhitnitsky [TO argues that 
the factorial growth (g 2 N c /2) 2k (2k)\(— of the large-order perturbative coefficients is 
related to instantons, thus being of non-perturbative nature. Important in this theory is 
dimensionality of the coupling constant. 

Gross and Periwal [^] (see item 12) proved that the /i-loop bosonic partition function 
(cut off in the infrared region) is bounded below and the perturbative expansion for the 
bosonic string diverges as Yl/h=o9 2h The series is not Borel summable, all its terms 
being positive. 

A look at the third column of Table 1 shows that most of the theories listed are 
characterized by an n\ large-order behaviour. This does not mean that all of them can be 
cured by the same resummation method: large-order behaviour is just one of aspects which 
determine the summation procedure. To each power series with coefficients listed in the 
3rd column of Table 1, there is a whole class of functions f(z) having the same asymptotic 
expansion. To specify the asymptotic expansion, one has to establish the angle (ray(s)) 
along which z approaches the origin; further, to pick out one function f(z) of this class, 
one has to add some additional information, according to the theory in question. 

As was already mentioned, these additional conditions are often formulated in terms 
of analyticity properties of the function expanded. 

In the next section, let us discuss the character of the additional conditions. 
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2 Asymptotic series. Summation methods 



If you are in doubt, 
examine the asymptotics 

How to deal with divergent series and under which conditions a power series can 
uniquely determine the expanded function are questions of fundamental importance in 
quantum theory. Power expansions are badly needed in physics, but additional input is 
required to ensure that they have precise meaning. These additional conditions should 
reflect some physical features of the system. 

We shall change the notation in this section, introducing the symbol z for the coupling 
parameter denoted previously by g. 

A formal series (^]) is called asymptotic to f(z) on the set S if the sequence Rn(z), 

N 

R N+ i(z) = f(z)-J2f n z n (10) 



satisfies the condition 

R N+1 (z) = o(z N ) (11) 

for all non- negative integers N and all z G S, where S is a subset of the complex plane 
having the origin as an accumulation point. 

The analogous notion of a series for which ( fLOD and (|TID are satisfied only for a limited 
number of N, N = 0, 1, 2, Nq, has also interesting applications. 

It is important that S is generally a subset of the neighbourhood of the origin. As 
an example to illustrate non-uniqueness, consider (^j for z approaching zero within the 
angle | argz| < (3 with (3 < tt/2, i.e., in the right-hand half of the complex plane. To a 
function f(z) satisfying ([!]), one can add a term of the form Aexp(— l/z a ) with A real 
and < a < tt/(2/3) without violating (|TT|) within S. 

A summation method states conditions under which an asymptotic series can deter- 
mine f(z) uniquely, and yields the explicit formula for f(z). 

2.1 Uniqueness by Borel summation 

The series 

oo 

E ^ (12) 
n=0 
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is called Borel summable if 

1) its Borel transform B(t), 

oo 

5(t) = $>XM (13) 

n=0 

converges inside some circle, | t |< S , 5 > 0; 

2) B(t) has the analytic continuation to an infinite strip of non- vanishing width bi- 
sected by the positive real semi- axis Re t > 0, and 

3) the integral 

1 roo 

g(z) = - e-V°B{t)dt, (14) 



z Jo 

called the Borel sum, converges for some z ^ 0. 

One can see the motivation of this summation method on a simple example. Consider 
a generic quantity D, calculated in perturbation theory with the coupling z, 

oo 

D(z) = Y,d n z n . (15) 

n=0 

This can be rewritten as 

°° POO 

D(z) = y d n z n (l/n\) / dte~H n . (16) 

If the series (|T3J) has a non- vanishing convergence radius p, the integration in (|TBp can be 
exchanged with the sum inside the circle. If we are outside the circle (or if the convergence 
radius is zero, p = 0) we exchange the order of integration and summation to define the 
series by the same expression (provided that the integral converges). In either case, we 
obtain 

D(z)= dte- t Yd n ^-= dte-'Bfzt), (17) 
Jo q n! Jo 

where B is the Borel transform of D. Taking d n = n\ as an example (finite-order coeffi- 
cients are irrelevant for the character of singularities) we obtain 

oo 1 

-t 1 flQ\ 



D(z) = / dte 



1 - zt 

This integral does not exist for z positive (where D(z) has a cut), nor is the Borel sum of 
such a series defined. For other values of z, the integral is convergent and selects one of 
the functions having Y^=on\z n as asymptotic expansion within the angle < argz < 2ir. 
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The summation can be denned in many ways; there are infinitely many functions with 
this asymptotic expansion. For instance, one can add any term of the form Aexp(— l/z a ) 
for < a < 1/2. 

Comparing the expansions ([12]) and ( |i~3l) we observe that the Borel transform has 
better convergence properties. To see this, let us compare the convergence radius p\ and 
p 2 of (|T2| ) and (|13[) respectively. We have 



1/pi = limsup y\a n \ (19) 



l/p 2 =hmsup W\a n \/n\. (20) 

n— >oo 

The factor \fn\ between p 1 and p 2 indicates that if the convergence radius of the original 
series is nonvanishing, that of B(t) will be infinite. Also, some singularities that are 
condensed at the origin will spread out to the complex plane. In a sense, the distribution 
of singularities of the Borel transform is a blow-up of that of the original function; the 
convergence properties of (p~3|) are better than those of (O). 



Nevanlinna [35] gave the following criterion of Borel summability (this criterion is a 
refinement of the Watson lemma ||3"6| , see |37]] and also for discussion): 

Let f(z) be analytic in JC(i]) defined by the inequality Re-j > - (with rj positive), a 
disc of radius \r\ bisected by the positive real semi-axis and tangent to the imaginary axis 
(see Fig. la), and let f(z) have the asymptotic expansion 

oo 

/(*) ~ £ a n z n . (21) 

ra=0 

If the remainder Rjsr{z) QlQl), after subtracting N — 1 terms from f(z), is bounded by the 
inequality 

\R N (z)\ < Aa N N\\z\ N (22) 
uniformly for all z G K{rj) and all above some value then f(z) can be represented 



by the integral ([14]) for any z G TC{rj). 

Taking QCD as an example, we see that the conditions of the Nevanlinna criterion are 
not satisfied, because the analyticity region is not a disc K.{rj) but a wedge of zero opening 
angle. 

2 This non-perturbative result (see ||] and 10 ) is obtained by combining asymptotic freedom with 
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Let us therefore consider some modifications of the Nevanlinna criterion. 



2.2 Analyticity vs. large-order behaviour: balance for unique- 
ness 

We shall briefly mention two modifications of the Borel method that are relevant for 
quantum theory. There must be a balance between large-order behaviour and A, the ana- 
lyticity and boundedness domain of f(z). The rule is simple: if the expansion coefficients 
grow faster than n\, one has to require a larger domain A to obtain a unique function 
from the expansion. If the large-order behaviour of the coefficients is tamer than n\, one 
can afford a smaller domain A to reach uniqueness in determining f{z). 

The Nevanlinna theorem can be generalized to the case that the opening angle of the 
analyticity domain is less than tt, in which case the disc IC(r]) transforms into a drop- 
shaped domain with its tip at the origin. 

Indeed, let F(z) be holomorphic in the "droplet" C(p) placed inside the sector | aigz\ < 
2^7, 7 < 1 with p < 7. We demand, in analogy with (|10D and (0), the bounds in the 
form 

N-l 

\f(z)-J2^z n \<Aa N (N\y\z\ N . (23) 

n=0 

The generalized Borel transform will have the form 

00 

B p (t) = J2a n t n /T(pn+l). (24) 

n=0 

We see that the condition (^) is more restrictive than (|22|) but, in compensation, its 
validity is required within a smaller angle. Correspondingly, the coefficients of the series 
B p {t) defined by (|4]) are less suppressed than those of B(t) defined by (|13"D (notice that 
B p {t) equals B(t) for p = 1). 

In discussing this generalized form for different values of 7 and comparing it with the 
original Nevanlinna theorem, we see that there is a balance between the size of the opening 
angle of A and the restrictiveness of the corresponding condition required for the series 
to determine the expanded function uniquely. This also means that Borel summability 
(i.e., p = 1, p(n) = n!)) plays no privileged role among the variety of possible summation 
methods. 

analyticity and unitarity of two-point Green functions in the complex momentum squared plane; see 
alsoQ. 
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In many practical problems, the Borel method nevertheless seems to be preferable, 
because most of the large-order estimates suggest an n\ behaviour of the perturbative 
coefficients (see Table 1). But this method simultaneously requires analyticity and the 
bound (^) in the z plane in an opening angle that is equal to tt. This is not always 
satisfied; in QED and QCD, as was already mentioned, the opening angle is zero in the 
coupling-constant complex plane, the origin being an accumulation point of singularities, 
and A being a wedge bisected by the positive real semiaxis and bounded above and below 
by circles that are tangent to it. 

This generalization of the Nevanlinna theorem gives no answer to such a situation; 
a zero opening angle would make us choose 7 and p equal to zero, but in this case the 
series ( p4|) for B p (z) would coincide with the original one ( |12"1) and there would be no 
resummation. Conditions for a unique determination of f(z) for A of such a shape have 
been obtained by A. Moroz f38| . They are given by the following theorem. 

Theorem 1 []38[| . Let f(z) 

1) be regular in the wedge W with the boundary |F(1/2;)| = F(l/R), F(z) = 
E~ * n M"). with 

POO 

H(n) = / exp(-e')t n dt, (25) 

Jo 

and continuous up to the boundary. 

Let the remainder R N (z), see flTU|), satisfy the bound 

\Rn(z)\ <Ah(N)\z\ n (26) 

uniformly in z G VV and N. 
Then 

00 

M(t) = £ U n /Kn) (27) 

n=0 

converges for \t\ < 1 and f(z) is uniquely represented by the absolutely convergent integral 

00 



f(x) = / exp(-e*)M(tx)dt (2* 







for any x G (0, R). 



The interested reader is referred to the original work [38] for details and also to a 



discussion in ||42j| . We just mention that, according to this theorem, the horn-shaped 
(zero angle) analyticity domain requires, for f(z) to be uniquely determined from the 
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f(z) at z — Uniform bound Transform Summation 

on Rn(z) 

1 Analytic £~ =0 a„^ n 

is convergent 

2 Singular, Ax^A!]^ B(t) = g(z) = 
opening angle = vr in /C(r/) and N > N £~ S£t n ± / °° e-'^B^dt 

3 Singular, Aa N (N\y\z\ N B p (t) = g p {z) = 
opening angle > in £(p), N > N £~ rggij JT e"* B p (t'*)dt 

4 Singular, A/x(JV)|*|* M(t) = g m (z) = M(tz) 
opening angle = in wedge W, N > N E^Lo ^n)^™ exp(-e*)dt 

Table 2: Analyticity vs. high orders: a balance is needed for uniqueness 

asymptotic series, a very tame large-order behaviour (given by (|26|) ) of the expansion 
coefficients, namely a at ~ (In N) instead of N\ required in the ordinary Borel summation 
method. Note the double exponential function in the integrand of (p5|), which causes this 
slow large-order behaviour. 

A survey of typical cases is presented here in the Table 2, together with conditions for 
a unique determination of f(z) from its asymptotic expansion. Conditions required for 
this are placed under one and the same item. Of course, one does not expect that f(z) 
will be determined uniquely from its perturbation expansion in real physical situations. 
Consequently, one cannot expect that the conditions placed in one item of Table 2 will 



be satisfied by a theory. This is indeed the case: QCD exhibits, according to 0] (see |10 
for general proof), analyticity in the wedge (see item 4 of the 1st column in Table 
2), but the large-order behaviour of the series is governed by the N\ rule (see item 2 of 
the 2nd column). This is no surprise: perturbation theory in QCD is incomplete and this 
discrepancy indicates the extent of non-uniqueness. The general rule for the use of Table 
2 is as follows: To ascertain uniqueness in the determination of a function f(z) out of its 
asymptotic series (|2l|) , one has to check if the conditions in all columns in one row are 



satisfied for the case considered; otherwise the function is either not uniquely determined 



'slightly modified by logarithmic terms 
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(b; 



(c) 





Figure 1: Three summation methods (see Table 2) for three different analyticity and 
boundedness domains: (a) the disc /C(r/), (b) the drop and (c) the wedge W. Crucial 
is not the size of the domain, but the opening angle at the origin. 

or may not exist .Q 

The disc JC(rj), the "drop" £(p) and the wedge W are depicted in Figs. 1 a, b, and c 
respectively 

Table 2 also shows that the large-order behaviour of the coefficients a n is not the only 
criterion of the (Borel or some other) summability of the series (pT|) . Even a series with 
a very tame behaviour of the coefficients a n may not be Borel summable, in spite of the 
radical suppression of the a n by the Borel factors l/n\. An example is discussed in [ j3£fl . 



2.3 Useful facts on power series 

Certain important facts on power series are sometimes overlooked in physical considera- 
tions. It may therefore be useful to recall some of them here because spontaneous intuition 
is often misleading. I apologize for some repetition of the text. 

1. The divergence of a perturbation expansion does not signal an inconsistency or 
ambiguity in the theory. An example to illustrate this is the result of Ogievetsky |J, who 
used the Borel summation method to sum a divergent perturbation series and obtained 
the exact result derived by Schwinger N5j, as is discussed here in subsection 1.1. 



4 The integral representations for g(z), g p (z) and g m (z) placed in the last column of Table 2 exhibit a 
symmetry which is revealed when appropriate substitutions of integral variables are made. 
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The problem is not that of convergence or divergence, but whether the expansion 
uniquely determines the function or not. The method of Feynman diagrams allows one to 
find, at least in principle, all coefficients of the perturbation series, which may determine 
the function uniquely even if the series is divergent, and may not do so even if it is 
convergent. This depends on additional conditions. 

2. The requirement fl2~T|) of asymptoticity of a perturbation series is not a formal 
assumption. It has physical content. 

When a perturbation series is divergent, one often re-interprets it automatically as 
an asymptotic series. This assumption is weaker than that of convergence, but is not 
self-evident. It is not only technical either. Its physical meaning is that there is a very 
smooth transition between the system with interaction and the system without it. 

For certain classes of observables the perturbation series is believed to be an asymptotic 
expansion. 

3. A very violent behaviour of the expansion coefficients a n at n —>■ oo might make us 
expect that no function with the property fl2~l~|) would exist. This fear is not justified; it 
was proved by Borel and Carleman (see for details) that there are analytic functions 
corresponding to any asymptotic series. 

4. Borel non-summability of a perturbation expansion alone does not signal an incon- 
sistency or ambiguity in the theory. Borel summability (or its lack) is not a fundamental 
classification criterion (see, e.g., |||). 

While large-order behaviour is often believed to be related to some very deep physics 
(see ||46||), the Borel procedure is just one of many possible summation methods and need 
not be applicable always and everywhere. There are power series with mildly growing 
coefficients which are not Borel summable; the problem is to find a method (not necessarily 
the Borel one) which is appropriate for the case considered. 

Example of a non-Borel summable series (although composed of very tamely behaving 
coefficients) which, nevertheless, is summable by the method based on Theorem 1 is 



considered in |!9| and discussed here in the end of subsection 4.1, (|56|). 

5. If the a n behave very violently at n — > oo (so that the Borel series J2^Lo a n z n /n\ has 
zero convergence radius) one might expect that it would be sufficient to replace n\ in the 
denominator by a sequence b n that grows faster than n\, in order to reach a more efficient 
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suppression of the a n . But this is done at a price: stronger conditions on analyticity are 
required for the summation procedure to be unambiguous. 

6. Asymptoticity is not a property of the series alone, but is contingent also upon the 
properties of the function expanded. These properties (e.g., analyticity) must be examined 
simultaneously with the asymptotic series, otherwise the same series can be summed to 
different functions. 

7. An asymptotic series need not be divergent. A wide-spread prejudice that gleams, 
often implicitly, in a number of papers is that an asymptotic series is always a divergent 
one. To see a quick counterexample, consider the function F(z) = g(z) + Ae~ a ^ z with A 
real and a positive, where g(z) is analytic at the origin. This is a two-parametric class 
of functions F(z) that all have, in the right half plane, the same asymptotic expansion, 
which is identical with the Taylor series for g(z); so the asymptotic expansion of all 
the singular functions F(z) is a convergent series because of analyticity of g{z). This 
convergent series is asymptotic to many singular functions, but its values converge only 
to one function, the analytic g(z). They do so inside the Taylor circle, which extends up 
to the nearest-to-origin singularity of g(z). 

3 Operator-product expansion and analyticity 
3.1 Operator-product expansion 

A standard way to investigate the interplay of perturbative and non-perturbative effects is 
the operator-product expansion (OPE) |^9|, [51]], in which the singularities of an operator 
product are expressed as a sum of non-singular operators with the coefficients being c- 
number functions. Being formulated before QCD, the operator-product expansion became 
a powerful (and actually the only) tool for examining non-perturbative effects in QCD. 
It consists in expanding the product of two local operators A(x) and B(y), for (x — y)^ 
tending to zero, into a series of local operators O n {Mx + y)), 

A(x)B(y) E C n (x - y)O n ((x + y)/2). (29) 

n 

where the C n are c-number functions. The meaning of the symbol ~ is discussed below. 
The operators O n are ordered according to their canonical dimensions, and their quantum 
numbers must match those of AB. As a consequence, the coefficients C n are ordered 
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according to their decreasing singularity degree, and only a finite number of them are 
singular functions infinite at the origin, while the remainder of the expansion is finite or 
tends to zero with vanishing distance. To any finite order in x — y, only a finite number 
of terms contribute. 

The short-distance behaviour of the C n is obtained, by dimensional counting, from 

C n {x) -> x d "- dA - dB (lnxm) p {l + 0(xm)) (30) 

with dAids and d n being the dimensions of A, B and O n respectively, and m being a 
quantity of the dimension of mass. In this way the canonical dimension of O n determines 
the singularity degree of the corresponding C n ; the operators O n with the smallest dimen- 
sions are dominant at short distances, being multiplied with the most singular coefficients 
C n . 

In momentum space, the expansion (|29"D takes the form (we put y — 0) 



i / dxe i<}X A(x)B(0) « C n (q)O n . (31) 
The symbol ~ is understood differently in different contexts. 

(i) While originally the validity of OPE was interpreted in the weak sense ||49|| , i.e., in 
that of sandwiching the product between two fixed states, from the very beginning it has 
been widely adopted ]5U) also to understand ( |2"5| ) and Q3"T| ) as (formal) operator relations. 

(ii) Different is the interpretation also in the perturbative and in the non-perturbative 
context. The standard derivation of OPE |49|, [51] relies on an analysis of Feynman graphs 
and remains within the framework of perturbation theory. The relations (^9[) and ([31]) are 
often understood as asymptotic relations (see, e.g., [0, ^,[5^]) which, in special situations 
(free- field theory, ||55||), may become convergent Taylor series. 



In application to non-perturbative phenomena, it is assumed that the left- and the 
right-hand sides of (|3T|) are approximately equal provided that a few first terms are taken. 
Then, at a critical dimension d r , the series breaks down and has to be supplemented with a 
term which is not of the Q~ 2k form. At higher values of k, other terms of different nature 
may appear. Relation ( |3~T]) can be understood as an approximate equality with finite 
number of terms on the right-hand side [|1^], with the hope that its validity and practical 
applicability (far from the point of expansion) can be established by phenomenological 
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analysis of experimental data. Application of this approach to hadron resonances has 
been very successful; agreement with experimental data and internal consistency of the 
QCD sum rules scheme are impressive. 

Non-perturbative effects manifest themselves also in the fact that, taking the vacuum- 
to- vacuum matrix element of fl3l|) , the terms on the right-hand side can be non- vanishing, 
whereas only the unit operator would survive within standard perturbation theory. 

(iii) Various attitudes are adopted to the topical problem of the validity region (angle, 
rays) of the expansion ([JI]) in the complex Q 2 plane. The expansion pT| ) reflects the 
properties of the left-hand side at the point of expansion and in a subset of its neigh- 
bourhood. In the euclidean (Q 2 > 0) region, the expansion is well defined in the sense 
that the separation of the large-distance contributions from the short- distance ones is well 
defined ||53|| . The Minkowski (Q 2 < 0) region, on the other hand, is most important for 
applications, because it contains the spectrum of the physical states. 

This indicates that the expansion ( pip can have different interpretations, according to 
the assumed physical content. These interpretations, however, have common mathemati- 
cal features, which we are going to discuss in the next subsection. We adopt the approach 
(which is mostly used in phenomenology and is formulated, e.g., in J56|) according to 
which the operator product expansion is originally built up in the euclidean region and 
then translated in the language of observables by analytic continuation to the Minkowski 
region. (For problems of the definition of euclidean Fermi and Bose fields, we refer the 
reader to [|57| and references therein). 

Note that along directions in the complex plane away from the euclidean ray different 
asymptotic expansions may apply or may not exist at all. 

Predictions in the Minkowski region are obtained by analytic continuation, which can 
be carried out under special assumptions of smoothness and stability. 

We shall now discuss conditions under which an asymptotic expansion can be extended 
to angles away from the euclidean semiaxis. It turns out that sign-definitness and inte- 
grability of the discontinuity along the cut in the Minkowski region play an important 
role. 
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3.2 Operator-product expansion away from euclidean region 

The combination of the operator-product expansion with analyticity represents a solid 
basis for the calculation of the QCD observables. Such calculations require knowledge of 
the properties of the expansion in the complex energy plane away from the euclidean ray, 
and analyticity allows one to establish them by means of analytic continuation. 

In practical applications of the operator-product expansion one typically meets the 
following generic situation. Let f(z) be holomorphic in C, the complex z plane (z = 1/Q 2 ), 
cut along the negative real semiaxis, z G [— oo,0]. Let f(z) be bounded by a constant, 

\f(z)\ <M (32) 

for z G C . Let the numbers a k , k = 0, 1,2, ...n — 1, exist such that the following n 
inequalities 

\R k (x)\<N k x k (33) 

are satisfied for all x = Kez, < x < d, with d being a positive number, n being a positive 
integer, and 

n-l 

Rn(x) = f(x)-J2a k x k . (34) 

k=0 

The problem is whether or in what sense the estimates (|33|) apply also to rays away from 
the positive real semiaxis. 

The condition ( |32j ) of boundedness by M can be generalized to cases in which the 
function of interest, F(z), is bounded by a polynomial in Q 2 = 1/z rather than by a 
constant; the theorem can then be applied either to F(z) multiplied by the corresponding 
power of z, or to the corresponding derivative of F(z). 

Notice that our considerations in this section, being equally well applicable to the case 
of a fixed, finite value of n as well as to the case of n infinite, are not restricted to the case 
of asymptotic series. We find the case of a finite n closer to physics because, first, one 
always knows only a finite number of the coefficients a k and, second, theorems dealing 
with a fixed value of n are more transparent and their application to the case of infinitely 
many inequalities (^) is straightforward. 

The solution to the problem very much depends on the character of the discontinuity 
of f(z) on the cut. 
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3.2.1 The case of integrable discontinuity 

1. Positive-definite discontinuity. Let us first assume that f(z) can be represented in the 



form of a generalized Stieltjes integral |59|j . 



/(*) = f tB-M, (35) 
Jo 1 + zt 

where the discontinuity is positive-definite all along the cut, i.e., p(t) is nonnegative for 
t > 0. Let p(t) be such that the integral in (^) exists for z complex, not negative, and 
that also the moments 

POO 

a k = / t k p(t)dt (36) 
Jo 

exist for all nonnegative integers k. 

The function f(z) has an asymptotic power expansion with the coefficients (— l) k cik, 

oo 

f(z) ~ V s (37) 

k=0 

for z approaching zero from the positive side, z —>■ 0+. Indeed, 



«-w = i->rj ifs rdi (38) 

and, consequently, 

|i?„(z)| =a n (^)|z| n <a n \z\ n (39) 

which, as expected, tends to zero one order faster than the highest power subtracted 

from f(z) (see (@)). We use the notation 

~a n (z) = r T^-r&t. (40) 
Jo \l + zt\ 

Analogous estimates are obtained for z approaching zero along a ray in the complex plane 
(but away from the Minkowski ray). In the right halfplane, certainly, a n (z) < a n , but 
this is not the case for Re z < 0, where the general bound on a n (z) is worse, due to the 
presence of the cut in this halfplane. The bound can be obtained as follows. We have, 
denoting z = re 1 ^, 

\l + zt\ = (l + 2rtcos0 + r 2 t 2 ) 1/2 . (41) 

Considered as a function of t at t > and <fi fixed, 1/|1 + zt\ has its maximum at 
t = to = — - cos0, where its value is 1/| sin0|. In this way we obtain ( ^9|) and 
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\Rn(z)\<—?—a n \z\ n (42) 



sm i 



for Kez > and Kez < respectively. Comparing (|42]) with fl39|), we see how the factor 
1/| sin 0| makes the estimate looser when the ray approaches the direction of the cut, i.e., 
when <p approaches zero. 

Possible applications to some QCD processes are discussed in subsection 5.2. 

2. Sign-indefinite discontinuity. The estimates are worse if the discontinuity along the 



cut is not positive definite. Physical relevance of such situations is pointed out in fl53"|. If 



the integral representation fl3"o| ) is still admitted, be it with p{t) sign-indefinite or even a 
complex-valued function on the integration interval, the resulting bounds fl39|) and ( |4*2"D 
have to be modified in the following sense. The function p(t) is assumed to be such that 
the integral (|35D is absolutely convergent and that the moments of \p(t)\ , 



oo 

k 



d k = / t k \p(t)\dt, (43) 
Jo 

are convergent integrals; this excludes wild oscillations of the discontinuity along the cut. 
The estimates (BUf) and fl4"2]) are replaced by 



\R n {z)\<a n \z\ n (44) 

and 

\R n {z)\<-^—a n \z\ n (45) 
| sm <p\ 

respectively. This may imply a considerable deterioration of the estimate relative to a 
sign-definite p(t): if the discontinuity on the cut oscillates, a n and |i? n (z)| will have very 
little to do with the expansion coefficient a n . 



Better estimates on \R n (z)\ than fl39|) , (f42|) and, respectively, fl44|) , (f45|) cannot be 



obtained unless some special assumptions are made about the explicit form of the function 
pit) in the respective cases. 

3.2.2 General case 

It is interesting to discuss a further generalization, although its relevance for physics may 
still be unclear. In the general case, when the integral ( p5| ) is not absolutely convergent or 
the moments ( |4"3"D are not defined, we can use a theorem by I. Vrkoc (see |5T]]) to obtain 
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an estimate. Assumed are holomorphy of f(z) and the bound (^) on f(z) to be valid in 
an angle 2a bisected by the real axis, and the inequality 

\R n (x)\ <N n x n (46) 

to hold for a positive integer n and for < x < d. When applied to a = tc (the values of 
n and d being fixed), the theorem implies that, along all rays between — 7r(l — 2~ l ) and 
7r(l — 2 l ), R m (z) satisfies the following inequality: 

\Rm(z)\ <Nm\z\ n/2 ' (47) 

for all | z | < do, where do is a positive number and I is a nonnegative integer. The 
coefficient N m in front of \z\ n ^ 2 ' is a linear combination of the |afc| for k running from 
m—[{n — l)/2'] + 1 to n — 1, [g] denoting the integer part of the real number q. 

This result is of relevance for expansions both with a finite and with infinite number 



of terms. Note that the bound (|46| ) on R n {x) for x positive implies, for z complex, the 
estimate ([|7]) on R m (z), where m is smaller than n. This resulting bound for complex z 
is considerably looser than that for x real, the number of exploitable terms being reduced 
from n down to m — [(n — 1)/2 Z ] + 1. This effect is the price paid for the extended 
validity region of (|47), and is becoming more and more pronounced with increasing I, 
i.e., with increasing angular deflection from the euclidean interval. We also see that the 
higher-order expansion coefficients (for k — m + 1, .., n — 1) cannot be used to improve 
the bound along a complex ray. 

The result shows quantitatively how quickly the expansion breaks down when the angle 
becomes large, i.e., when one approaches the cut on the Minkowski side of the real axis. As 
is usual in analytic extrapolations f6"0f , a point-by-point continuation from approximate 
data is impossible, and an extrapolation up to the cut can be safely done onto an interval 
but not to a point. The procedure should be based on statistical analysis of experimental 



data, pT 



Since the estimate fl47|) might seem rather loose (note that n in the exponent on the 
right hand side is divided by 2 l ), it is interesting to look for a function which would 
saturate it. A set of functions saturating (f47|) for different nonnegative integers n can be 
generated by using the function exp(—ia(ln z) 2 ) with a real.0 
5 I am grateful to I. Vrkoc for pointing out this to me. 
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Further modifications of the results of this section will arise when the logarithmic 
dependence of the expansion coefficients C n (x) indicated in (|30"D is taken into account. 
Physical relevance of this generalization is evident. 

4 Renormalons and instantons 
4.1 Singularities in the Borel plane 



The integral (13) defines the Borel sum of the series ( |T2"D if the conditions 1), 2) and 3) 
from subsection 2.1 are satisfied. 

The condition 1) would be violated if the a n were to grow faster than n\. As follows 
from Table 1, this is not the case in typical situations. 

We generally do not know the nature or distribution of singularities to assess the 
validity of the condition 2). Analyses of large-order behaviour indicate that 2) is violated 
by an infinite set of singularities located on the real axis (renormalons, instantons). Then, 
the integral (|T4] ) is not well-defined and should not be used. In quantum theory, however, 
it is nevertheless used, with the aim to reduce the problem of summation to that of 
defining the integral QT4] ) by means of a suitable modification of the integration contour. 

The spectrum of ambiguitites is wider, however. By introducing the renormalon lan- 
guage we consider only special sets of equidistant singularities (poles, branch points) 
whose parameters have physical interpretation. Those lying on the positive real semiaxis 
produce ambiguity and remind us that non-perturbative input has to be added. 

The electromagnetic current-current correlation function is a useful example to discuss 
a typical structure of singularities in the complex coupling-constant plane. Denoting this 
function IP", 

W v = i J d 4 xe" i9a; (0| T{f{x)f{0)) |0> (48) 
= (ff"V - gVM-g 2 ) (49) 

and taking R, the ratio of the total cross section for e + e~ — > hadrons to that for e + e~ — > 
which is related to the imaginary part of II, 

R(s) = 127rlm n(s + i0 + ) , (50) 
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we introduce a modified quantity D(Q 2 ) defined as 

D (Q 2 ) = -^ 2 g 2 H(Q 2 ) (51) 

(where Q 2 = —q 2 , i.e., Q 2 > in the euclidean domain), to avoid inessential logarith- 
mic terms and the dependence on the renormalization scale \i ||64|| ). The perturbation 
expansion of D in powers of the coupling constant has the form 

D ^ 1 + ^lf^ Dn{as{Q 2 )r _ (52) 

* n=0 

The dependence of D on a s (Q 2 ) exhibits a complex structure of singularities in the 
coupling constant plane at and around the origin |§. Their nature can be conveniently 
displayed by studying the corresponding Borel transform. The behaviour of the running 
coupling constant g(fi 2 ) is determined by the j3 function, (3{g 2 {ji 2 )) = yU 2 (d/d/z 2 )g 2 (/x 2 ), 
which has the perturbative expansion (3{g 2 ) ~ —bog 4 — big 6 — ... . 

The structure of the singularities of the Borel transform that are nearest to the origin 
can be seen from the following (simplified, see |62|| ) expression for D 



Q 2 J dk 



2 k 2 a s (k 2 ) 

(P + Q 2 ) 3 ' 1 ' 

which at the infrared and the ultraviolet end of the spectrum has the following behaviour 

a s {Q 2 )Y.[ba s {Q 2 )/{^)fn\/2 n+1 (54) 

n 

and 

a s (Q 2 )T,[-^s(Q 2 )/W} n n\ (55) 

n 

respectively, where b is the first coefficient in the Gell- Mann-Low function. This corre- 
sponds to an n! large-order behaviour, with the nearest singularities located in the Borel 
plane at 8n/b and —Amjb respectively. 

Perturbation theory suggests the following structure of the singularities of the Borel 



transform [Bl 



(1) Instanton-anti-instanton pairs |fl8| , |I~9|j generate equidistant singularities along the 



positive real axis starting at t = 4, for t = 41,1 = 1,2,.... Balitsky |J5] calculated the 
behaviour of i? e + e -^hadrons near t — 4 and found the leading I — I singularity to be a 



28 



branch point with the power ^(iVf — N), where N and Nf is the number of colours and 
of flavours respectively. 

(2) Ultraviolet renormalons are generated by contributions behaving as c^+i ~ (—bo/l) k k\ 
for / = 1,2,..., leading to singularities located at t — —2l/b on the negative real axis 
(where b = 87r 2 &o), with a branch point for 1 = 1. 

(3) Infrared renormalons are generated by contributions behaving as c^+i ~ (b /l) k k\, 
I = 2,3, leading to singularities located at t = 2l/b. Near the first of the points, the 
singularity behaves [|ll|] as (bt — ^i) -1-2 *' 60 , \ = bi/bo. By absence of a dimension-two 



condensate the existence of the I = 1-singularity is excluded. 

Other singularities are not assumed to exist in the Borel plane. But there may be 
still another source of ambiguities, due to possible singularities at infinity. They may 
be non-Borel summable, no matter how tame the large-order growth may be! A trivial 
example is the series S^Lo zn ■> whose coefficients are constant. It has no singularities in 
the Borel plane, but the series is Borel summable only for Rez < 1, because otherwise the 
Borel integral does not exist. A more sophisticated example is Y^=oC n z n with 

c n = c n fi(n) 

roc 

H{n) = / exp(-e') t n dt (56) 
Jo 

c = \ c \e w , < < tt/2 

(see ref. |39|, Example 3). Note the double exponential in the integrand, which makes 
and c n rise extremely slowly. At first sight, this series should be Borel summable, even 
more that there are no singularities in the Borel plane, the convergence radius (condition 



2)) being infinity. But the condition 3) is not satisfied: the integral {TA) does not exist. 
Indeed, 



dte-'/'E^" (57) 

J u n=0 "'■ 

is divergent, because B(t) grows faster than e* with t — ► +oo. But the series Y^=o c nZ n is 



summable by Moroz's method |38 



4.2 A further generalization of Borel transformation 

The functions B p (t) and M(t) defined in Table 2 are generalizations of the Borel transform, 
which can be used in the various situations listed in Table 1 to reduce non-uniqueness, 
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provided some additional information is available. More about the properties of B p (t) and 



M{t) can be found in 



42| and in references therein. 



I will now discuss another type of generalization of the notion of Borel transform, p3 



pE7|| , which makes use of the specific structures mentioned in the previous subsection. 
Let me first make a general remark. After having exposed mathematical methods in 
sections 1 and 2, we have passed to practical aspects of the operator-product expansion 
(estimate of the remainder after the n-th term, see section 3), and of the assumption (in 
section 4) that two-point Green's functions have special singularities in the Borel plane, 
the instantons and the renormalons, a structure that is now almost universally adopted. 
Existence of these structures in the Borel plane is based on various mathematical models 
developed in late 70's and early 80's, but nowadays they are mostly considered as true 
features of Nature. 

Brown, Yaffe and Zhai and Beneke |)3, 0, |47j use information about the structure of 
the first infrared renormalon. Assuming first the special case of a one-term (3 function, 

P{g 2 ) = -b Q g\ (58) 

expanding Im II and II in powers of the coupling constant with the expansion coefficients 
a n and c n respectively and defining their respective Borel transforms 

na r 



Ait) = £ 



jr(n + i; 



and 



C(t) =5 + E 



nc r , 



n=iT(n+l) 



(59) 



(60) 



Brown and Yaffe obtain, by comparing the expansion coefficients, the following rela- 
tion between A and C: 

A(t) = sin(b t)C(t) , (61) 



which turns out to be a consequence of renormalization-group invariance |47| . Defining a 
modified Borel transform ^(t) by 



(62) 



^ T(n + l + Xty 

(thereby accounting for the first infrared renormalon), the authors of |4^J consider 
the case of a general beta function, chosen in such a scheme that its inverse contains two 
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terms: 

l/(3(g 2 ) = -l/(bog 4 )+A/(b g 2 ). (63) 

They find that, for this form of f3(g 2 ), the relation remains valid also for A(t) and 
C(t), the modified (according to (|62|)) Borel transforms of Imll and II respectively: 

A(t) = sin(6 t) C{t) . (64) 



In this way, the endeavour to retain the relation (^H]) also for the case of a general f3 
function motivates the authors to introduce the generalized Borel transform in the form 

flU). 

It was already pointed out that the concept of renormalon is at present applied to con- 
crete physical situations. This poses a topical problem: To what extent are renormalons 
physical concepts and to what extent are they just artefacts of our, maybe inadequate, for- 
malism? Some authors argue that the series of renormalon- type graphs is ill-defined, 



or that infrared renormalons might be an attribute of the renormalization scheme used. 
The Borel plane formalism is, nevertheless, successful in phenomenology. 

4.3 Phenomenology in the Borel plane 



Soper and Surguladze |67J use information on the singularities in the Borel plane to 



compute the ratio R of the width for Z— > hadrons to that for e + e~. The result is 
impressive. Considering R as a function R(s) of s, the cm. energy squared of the e + e~ 
annihilation, one obtains the measured R as i?(M|). Representing R(s) and D(s) as 
Ro[l + 7l(s)} and D [l + T>(s)] respectively with 1Z and V expanded in powers of the 
coupling constant a s (s)/n, the authors of |67j use, in view of the relations (|50|) and (pl|), 
the formula 

K(s) = r d9V(se w ) (65) 

to calculate 7Z{s) from T>(z). Both TZ and D are expanded in perturbation series to third 
order, where the coefficients are explicitly known. Assuming now a standard large-order 
behaviour of the perturbative expansion of V, the authors write the Borel representation 
for it, 

roo 

V(Q 2 )= dzex P (-nz/(a s (Q 2 ))V(z), (66) 
Jo 
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where the Borel-transform coefficients T> n have better large-order behaviour than the 
D n , those of the original series. Combining (|65| ) with fl6"6|), one expresses 1Z directly in 
terms of T>(z) through an integral over z from to oo, in which the upper limit can be 
safely replaced by z m3iX , with (3oZ max around 0.5 or higher. The ultraviolet renormalon at 
(3 z = —1, being the singularity closest to the origin, controls the large-order behaviour 
of the perturbative series; the authors of |67| move it farther by means of the conformal 
mapping 

y/l + p Z + 1 

thereby improving the convergence rate of the power series in the Borel plane. Note 
that the other set of dangerous singularities, the infrared renormalons (the nearest of 
which is at f3$z = +2), are not pushed away from the origin by (|6TD, but on the contrary 
become closer to it. Although the infrared renormalons could be pushed away by means 
of the optimal conformal mapping [68] (which leads to the fastest convergence), Soper 



and Surguladze |57| prefer to soften the nearest infrared renormalon by multiplying D 
with its power, because they know the numerical value of the exponent. Making good use 
of this valuable piece of information, they reach a high accuracy in determining a s (M|), 
estimating that the uncertainty arising from QCD is less than 0.5 per cent. 

4.4 Renormalons and OPE 

A scheme to use renormalons to construct non-perturbative terms in the operator-product 
expansion has been proposed by Grunberg p9[ . Defining from H(Q 2 ) the renormalization- 



group invariant quantity 

F?(0 2 \ = 

dlnQ 2 dlnQ 2 



R(Q 2 ) = -r^--^Uo, (68) 



he uses the OPE representation 

R{Q 2 ) = R PT {a) + G {a) + (69) 

where a = a(/j,) is the coupling at scale /x, Rpx(a) = Y^=o r n(y n+1 the perturbative contri- 
bution and Go(a) is the leading condensate contribution. For Go (a), the renormalization 
group equation yields 

G {a) = C{fi 2 /Q 2 ) d/2 exp{d/{2(3 a))a 5 [l + 0{a)}, (70) 
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where 5 is expressed through the /3-function coefficients, and the /^-dependence of the 
factor (/i 2 / Q 2 ) d / 2 is compensated by that of the subsequent exponential function. 

The function Rp^{a) is then approximated by first infrared renormalon singularity. 
Using the Borel representation 

POO 

R PT ( a )= dzexp(-z/a)Rp T (z) (71) 
Jo 

with Rp T (z) approximated by the first renormalon Rp T0 (z), 

R*(z) = K(fi 2 /Q 2 )- Z ^(l - z/z )^[l + 0(1 - z/zo], (72) 



where K is a scale-independent factor, Grunberg arrives [^9j at 

lmR PTfi (a) = ±^ 7 r(l - 7 ) sin(7r(l - 7 ))(/i 2 / 'Q 2 y z ^ exp(-^ /«)« 1 - 7 [l + 0(a)}. (73) 

The sign depends on the choice of the integration contour; thus, the resummed Rp^(a) 
cannot determine R uniquely. Consistency requires the ambiguity to be removed by some 
non-perturbative input. Assuming C in (f70[) complex, Grunberg |69| compares (|73|) with 
dTO|) to find relations connecting zq/3q and 7 with d and 5 and to fix the position and the 
power of the singularity. 

This method allows non-perturbative terms in the operator-product expansion to be 
constructed from perturbation theory. The author gives a convergent scheme to deal with 
perturbation theory in presence of infrared renormalons. The procedure can be extended 
to non-leading IR renormalons. 

Grunberg's result shows that there are "non-perturbative" terms in the operator- 
product expansion which have their origin in perturbation theory. The author considers 
also the possibility that all "non-perturbative terms" in the expansion, taken to all orders, 
could be determined in principle from perturbation theory. This interesting possibility 
could be a good starting point for approximations, although it does not mean that all 
QCD could be fully described in terms of perturbative concepts. Within perturbation 
theory one can account for the infrared domain up to terms of the form 

exp(m/(/? a(Q)) ~ (A 2 /Q 2 ) m (74) 

(see a discussion in [0), where m is the order of IR renormalon. On the other hand, 
an additional non-perturbative input may also have implications in the Borel plane. The 
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presence of renormalons can signal deficiency of the perturbative approach and indicate 
a nonperturbative supplement. 



5 Further applications 

5.1 Resummation of bubble chains 

The renewed interest in calculating higher-order perturbative corrections and in examining 
the high-order behaviour of perturbative series is intimately related to the investigation of 
renormalization scale and scheme dependence of a truncated series as well as to attempts 
to estimate its uncalculated remainder. In the past various criteria for finding a suitable 
renormalization prescription were proposed; they are based on an estimate of the size of 
the remainder. Examples are Stevenson's principle of minimal sensitivity ||, Grunberg's 
notion of effective charge [7I|, and the BLM method of scale setting [[72] by Brodsky, 



Lepage and Mackenzie. The BLM prescription, which is based on an analogy with QED, 
was further developed recently. It is a method of estimating higher-order perturbative 
corrections of a physical quantity, provided that the first approximation is known. It 
consists in the use of some "average virtuality" as scale in the running coupling. Instead 
of working with fixed scale, 

a s (Q 2 ) [ d 4 kF(k,Q), (75) 



one averages over the logarithm of the gluon momentum: 

<QIlm) J d'kF(k, Q) = a s (Q 2 ) J d*k - ^a s (Q 2 ) In F{k, Q) . (76) 

This replacement amounts to accounting for higher-order terms in powers of a s (Q 2 ) by 
making use of the renormalization-group evolution 

a s (-k 2 ) = a s (Q 2 )f(^^) n ~\-H-k 2 /Q 2 ))^ (77) 



n=l 



(in the leading-logarithm approximation) and retaining only the first two terms in the 
sum. This approach was generalized |7^, [80], |ST] by introducing the running coupling con- 
stant a s (k 2 ) directly into the vertices of Feynman diagrams, with k being the momentum 
"flowing" through the line of the virtual gluon. This modification means replacement of 
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by 

u s (Q* 2 ) J d 4 kF{k,Q) = J d*ka s (-k 2 )F(k,Q) (78) 

with a s (x 2 ) = An / (P Q \n(x 2 / Aq CD ). Note, however, that the beta function f3(a s {Q 2 )) is 
approximated by its first term only. 

This method has been applied in phenomenology to various QCD observables, such as 
r decay hadronic width and heavy-quark pole mass [^, ^T[, semileptonic B-meson decay 
|82f and the Drell-Yan process [jS3|l . It makes maximal use of the information contained in 
the one-loop perturbative corrections combined with the one-loop running of the effective 
coupling, thereby providing a natural extension of the BLM scale-fixing prescription. 

Methods generalizing the scale-setting procedure developed by Brodsky, Lepage and 
Mackenzie met with criticism H, RH| showing that the summation based on "naive non- 



abelization" is renormalization-prescription dependent. Further research will clarify the 
issue; some of its aspects are discussed in Sec. 6. 

A resummation that is renormalization-scheme invariant in the large- Nf approxima- 
tion has recently been proposed by Maxwell and Tonge ||86|| . It is based on approximating 
the effective charge beta-function coefficients by the part produced by the highest power 
of b = (11N C — 2Nf)/Q, and accounts for exact all-orders results in this approximation. 
This method allows one to include, in any renormalization scheme, the exact perturbative 
coefficients not only at NLO but also at NNLO. The approach is used to critically examine 
the reliability of calculations of fixed-order perturbation theory for various observables and 
sum rules. Application of the RS-invariant resummation to different physical processes 
indicates varying reliability of the standard fixed-order perturbation theory. 

Ellis et al. |87| use Pade approximants to develop another method of resumming the 
QCD perturbative series. The authors test their method on various known QCD results 
and find that it works very well; in particular, they reach a weaker renormalization- 
scale dependence of the calculations. The problem why Pade approximants reduce the 
renormalization-scale dependence has recently been considered in ref. [J88[ . The author 
finds that in the "large (5q limit" (when the beta function is dominated by the 1-loop contri- 
bution) diagonal Pade approximants become renormalization-scale independent, thereby 
sharing an important feature of the exact result. 
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5.2 Use of operator-product expansion 

Operator-product expansion has proved to be very useful in practical calculations in con- 
nection with the tau decay hadronic width, the heavy quark pole mass, heavy-light quark 
systems, semileptonic B-meson decay, heavy quarkonia and the Drell-Yan process. Re- 
newed interest in the evaluation of power corrections to QCD predictions has resulted 
in attempts to understand the implications of the presence of renormalon singularities in 
these processes. 

Starting from analyticity and operator product expansion Braaten, Narison and Pich 



39| used the measurements of the r decay rates to determine the QCD running coupling 
constant at the scale of the r mass M r =1.784 GeV. The ratio R T , 

T[t- -> z/ r hadrons(7)] 
T[t -> z/ T e i/ e (7)J 

(where (7) represents possible additional photons or lepton pairs), is represented in the 
form 

r M r Hs 
127T J (1 - s/M 2 ) 2 (l + 2,/M r 2 )Imn( S )^ (80) 



where n(s) is a combination of the two-point correlation functions for the vector V£j 
ipjl^ipi and axial vector Afj = ipj'j^'j^ipi colour singlet quark currents, with coefficients 
given by the elements V u d and V us of the Kobayashi-Maskawa matrix, the subscripts 
i, j =u,d,s denoting light quark flavours. 

This integral cannot at present be calculated from QCD, because the hadronic func- 
tions are sensitive to the non-perturbative effects that confine quarks in hadrons. But one 
can make use of the properties of the correlating functions, which are known to be analytic 
in the complex s-plane cut along the positive real semiaxis. This allows the integral in 
fl8"0D to be expressed as a contour integral along the circle of radius M 2 : 

67ri£ |=M2 (l - s/M 2 ) 2 (l + 2 S /M 2 )n(,)^ . (81) 

Representing now Tl(s) as the operator product expansion over local gauge invariant 



operators, the authors of []89| point out that the double zero in the kinematic factor of 
the integrand at s = M 2 suppresses the contribution from the dangerous segment close to 
the branch cut, where OPE has little chance to represent appropriately the function. (As 
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a further means of suppressing bad influence of this segment, integral averaging smeared 
over an energy of the order Aqcd can be used.) Application of the estimates discussed in 
subsection 3.2 will give quantitative background to this qualitative assumption, specifying 
how reliable the expansion is at the points \s\ — M% that are far both from the euclidean 
interval and from the immediate neighbourhood of the point s = M^, thereby possibly 
modifying the accuracy of some of the estimates. 

The estimate of a s and of the QCD condensates from e + e~ — > I — 1 data was re- 



examined in |5D] using r-like inclusive process and QCD sum rules. The resulting values 
confirm previous sum rules estimate based on stability criteria. Implication on the value 
of ot s from r-decays gives 0.33 ± 0.03 as the value of a s (M T ). 

Applications of the operator-product expansion to heavy quark physics (such as B and 
Ab decays, the problem of heavy-quark pole mass) are based on an extension of the idea of 
the operator-product expansion to a new situation, when the expansion is made in terms 
of operators in the heavy-quark effective field theory, in inverse powers of a quantity of 
the scale of the heavy quark mass tuq [ 7(| . Heavy-quark field theory is, to leading order, 
a theory for light quarks in the field of a static colour source. Evidently, mQ, this key 
quantity, is not uniquely defined because quarks are unobservable; the so-called pole mass 
r?7,Q 0le (i.e., the position of the pole in the quark propagator) is well-defined in each finite- 
order perturbation theory |S| but, as the authors show, no precise definition of the pole 
mass exists once nonperturbative effects are included. They argue that any consistent 
definition of this quantity suffers from an intrinsic uncertainty of order [pole mass times 
Aqcd/^q]- This circumstance can be described in terms of an infrared renormalon of the 
type ([54]), which produces a factorial divergence of the perturbation series. 

As in the r decay case, also in the case of the heavy-quark systems the approach 
starts from the operator-product expansion in euclidean region, this time in the kinematic 
invariant v.q, where v is the velocity of the heavy quark. Again, the problem leads to a 
contour integral of the type (|81|), which has a segment close to the cut, where the expansion 
is not safe, and it is expected that smearing over an energy of the order Aqcd would make 
valid the OPE computation of the differential distributions. Results for semileptonic B 
and A b decays are obtained in |7I| . 

The relation of the heavy quark mass to the expansion parameter is investigated in 
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|70| . It is shown that the heavy-quark effective theory suffers from a series of ultraviolet 
renormalons which are not Borel-summable. The heavy-quark effective theory should be 
independent of the heavy-quark mass; since the S'-matrix elements in a confining theory 
have no poles corresponding to a physical quark there is no natural choice of the expansion 
parameter. The ambiguity of perturbation series is related to an infrared renormalon in 
the pole mass and expresses the necessity of including an additional (and ambiguous) 
mass term. 

The structure of renormalons in the heavy-quark effective theory is examined in | |77| . 
The authors show that the way in which renormalons appear in effective theories is not 
universal, but depends on the regularization scheme used to define the effective theory. 
In the case of kinetic energy operator, they find that the leading ultraviolet renormalon 
is absent in all but one regularization schemes. 

The authors of develop a nonperturbative method for defining the higher-dimen- 
sional operators of the heavy quark effective theory such that the matrix elements are free 
of ambiguities due to the UV renormalon singularities. They define a "subtracted pole 
mass" from which the renormalon singularities are subtracted in a nonperturbative way 
and whose inverse can be used as the expansion parameter. 



For a critical review of recent achievements in the heavy quark theory see ref. [79 



6 Renormalization ambiguities 

Results of any finite-order calculation depend on the renormalization prescription, al- 
though the theory as a whole must be invariant with respect to the renormalization 
group. This issue has been lively discussed in QCD since many years in connection with 
theoretical predictions of observable quantities. As the renormalizat ion-group dependence 
of a finite-order result can never be fully removed, the effort has been focused on looking 
for such a renormalization prescription (RP) and such a method of approximation which 
would possibly minimize the dependence on the choice of the RP. The extent of ambigu- 
ities connected with renormalization-group invariance is often not fully appreciated; we 
shall therefore explain the main concepts and the actual origin of the ambiguities. In 
expounding the subject we shall follow the way outlined in ref. f83}| ; for practical aspects 
of the renormalization scheme dependence of the perturbative QCD approximations see 
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also [0, 0. 

Let r(Q) be a generic observable quantity that depends, for simplicity, on a single 
external momentum Q, 

r(Q) = a [l + ri a + r 2 a 2 + ...] . (82) 

On the right-hand side, the Q-dependence of r(Q) is hidden in a Q dependence of the 
expansion coefficients r k . Besides this, both a (the renormalized coupling constant, often 
called coupling parameter or couplant in this context, to emphasize its running character) 
and the coefficients r k depend on the renormalization prescription (RP), which consists 
of the choice of the renormalization scale fi (of dimension energy) and of the so-called 
renormalization scheme (RS). Thus, 

a = a(/i, RS) 

r k = r k (Q/fj,,RS). (83) 

The quantity r(Q), on the other hand, should not (by definition) depend on fi or RS, 
because it is an observable quantity^ 

One can distinguish the following stages of fixing the renormalization-group ambigui- 
ties. 

1. Renormalization scale. The dependence of the couplant a on the scale fi is governed 
by the differential equation 

da [^ RS) = /3(a) = -ba\fi, RS) (1 + caU, RS) + c 2 a(//, RS) + ...). (84) 
din ix 

The choice of the renormalization scale consists in selecting, when equation fl84|) has been 
solved, the value of the parameter /i, with respect to which the theory is invariant. 

2. Renormalization convention. While b and c, the first two coefficients on the right 
hand side of ([84]), are uniquely determined as functions of Nj and iV c , 

b = (llN e - 2N f )/6 
c = (51iV c - 19N f )/(llN c - 4N f ) 



6 Strictly speaking, this statement is true in the fictitious world of perturbative physics. Genuinely 
independent is the sum r(Q) + r NPt (Q), where the second term is the nonperturbative part. But such a 
statement has no predictive power unless something is known about r NPt (Q). The standard approach is 
to assume that r(Q) and r NPt (Q) are separately independent of the renormalization prescription. 
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Symbol Name to be fixed what happens 

R-scaie renormalization scale \x scale chosen 

RC renormalization convention Cj form of curves chosen 

RS renormalization scheme Cj, Ars curve chosen 

RP renormalization prescription Cj, Ars, A* complete choice 

Table 3: Stages of renormalization fixing: renormalization is fully fixed by specifying 
the renormalization prescription 

(where Nf and N c is the number of massless quark species and the number of colours 
respectively), all the other coefficients q are completely arbitrary and define the renor- 
malization convention, see Table 3. One can consider different renormalization conventions 
(i.e., different choices of the q), among which 't Hooft's convention ||, in which all the Cj 
vanish, is remarkable for its simplicity and is suitable for general considerations; in this 
convention, the right-hand side of fl84]) consists of two terms, while (|32D is generally an 
infinite series. One can think of another special choice, in which fl8"2|) is a truncated ex- 
pression and ( p4|) is an infinite series (Grunberg's convention, effective charges approach, 

TO- 

3. Choice of the renormalization scheme. Once the coefficients Cj are chosen and some 
boundary condition on a(fi, RS) is specified, the equation (|4|) can be solved. Concerning 
the different ways of specifying the boundary conditions on a(/x, RS), a popular way con- 
sists in selecting the scale parameter Ars, defined by the following implicit equation [|J 
for the solution of the differential equation fl84|) 



A 4 1 . ca r a , 
b In - — = — hem h / ax 

Ars a 1 + ca Jo 



1 

+ 



\BW)(x) x 2 (1 + cx) 



15) 



where B^ n '{x) = 1 + cx + C2X 2 + ...c n _ix n_1 . The integrand on the right-hand side of ([85]) 
is regular at the lower integration limit and vanishes in 't Hooft's convention, in which all 
the Ci vanish. 

According to (|S~5{ ), the choice of Ars amounts to selecting a special curve among the 
solutions to (0), which in 't Hooft's convention has a particularly simple interpretation: 
Ars is that value of fi for which a(/x, RS) becomes infinite. Equations similar to ([85]) 
govern also the dependence of the couplant a(/i, q) on the parameters q; the equations 
contain the same fundamental (3 function and introduce no additional ambiguity. 
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To sum up, the full size of ambiguities connected with renormalization group imparl- 
ance involves the choice of the renormalization scale /i, the choice of the values of the 
Q-coefficients, and the choice of the boundary condition on a(p, RS) by specifying the 
renormalization scheme. A survey of degrees of renormalization fixing is given in Table 3; 
let me just point out that the terminology (convention, scheme, prescription...), although 
widely adopted at present, is sometimes used in a different way. 

The explicit dependence of the expansion coefficients on the renormalization pa- 
rameters /i, Ars and the Cj is given by the condition that r^ N \Q/p, RS), the sum of the 
first N terms of the right-hand side of (|82|), 

2V-1 

r W (Q//i,RS) = r k {Q/p,,R$)a k+1 {p,R$), (86) 

k=0 

when differentiated once with respect to In /i, q and In A RS respectively, should be of the 
order 0(a N+1 ) (note that, in the latter case, it is not the value of Ars what varies at fixed 
RS, but the change consists in a transition from one renormalization scheme to another). 
This allows us to introduce the renormalizat ion-group invariants pi, i = 1,2,3,... and 
represent the coefficients in terms of them in the following form 

niQ/^RS) = 6 In/- -MQ/Ars) 

r 2 (Q/fM, RS) = p 2 - c 2 + r\ + cr x , (87) 

etc., where only p\ depends on Q, while the pi, % = 2,3, ... are numbers, the dependence 
of p\ on Q being obtained from fl87D by setting p = Qtf\ 

Thus, although the expansion coefficients r« and q are RP-dependent, certain combi- 
nations of them, the function p\ and the numbers RP-invariants. Their 
knowledge is relevant for the the extraction of QCD observables from experiment (see, 
e.g 



84] , [91] , Pq , |96|, |97[). There are two different sources of uncertainties in the theoretical 
determination of a QCD quantity, one being connected with higher-order (uncalculated) 
perturbative corrections, the other with the renormalization scheme and scale ambiguity. 
Since the exact expression (perturbative sum) must be RG invariant (see the footnote 



It might seem paradoxical that the function pi, while being a renormalization-group invariant, is 
represented here as a function of Q/Ars, where Ars is n0 invariant. As p\ is dimensionless, its Q- 
dependence takes place only through dimensionless quantities. Then, to preserve invariance, the symbol 
pi must represent different function prescriptions for different schemes. 
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6), there is a widespread belief that a weak RG dependence is a signal of a small large- 
order correction, i.e., that the strength of the RG dependence of a particular observable 
is related to the size of its higher-order corrections which, by this, can be minimized by 
choosing the renormalization prescription that has the weakest dependence on the RG 
parameters. 

7 Concluding remarks 

A typical feature of the present status of the QCD perturbative corrections is the trend to 
avoid explicit calculation of higher-order corrections and, instead, to improve the result 
by making use of some additional information. Such information may be, for instance, 
the renormalization group invariance (which allows us to introduce the running coupling 
constant instead of the fixed one), analyticity, or the structure of singularities in the Borel 
plane. 

As was already mentioned, the notion of renormalon, originally introduced and used 
to investigate interesting mathematical models, is now widely believed to have concrete 
background in physical phenomena. This wide-spread opinion meets also with scepticism, 
which argues that some singularities in the Borel plane might be an attribute of the 
renormalization scheme used p6| . 

It seems that the present effort in further developing the idea of renormalon and the 
corresponding formalism will be helpful in finding a language appropriate for supplying 
the non-perturbative sector with new physical ideas. Generalizations of the scale-setting 
procedures become valuable by implementing new physical information without the need 
to calculate higher-order perturbative corrections (which is not only cumbersome but 
also doubtful, if the series is divergent or convergent with a bad rate), by using some 
additional, perturbatively independent, information. This idea is not new, appearing in 
theoretical physics whenever technical difficulties force one to look for methods allowing 
the exploitation of all information on the system, including whatever is not adequately 
taken into account by the existing formalism. 

This additional information is needed not only whenever one is unable to calculate 
higher-order corrections. It seems evident that perturbative quantum chromodynamics is 
not self-sufficient on principle, even if we had calculated all the perturbative coefficients 



42 



of any order. The asymptotic perturbative series calls for new dynamical input. The 
operator product expansion, however bringing valuable physical contribution, is not suffi- 
cient either, being as a power series connected with a similar ambiguity. It is intriguing to 
have the missing input supplied by the general principles of the S'-matrix theory, among 
which analyticity occupies the central position. This is not only due to its rich physical 
content, but also because, as was pointed out in subsection 1.2, the S'-matrix possesses 
the standard analyticity properties even for gauge theories with confinement []nj. This 
fact distinguishes all summation methods that are able to make use of analyticity, as 
additional input, to remove the ambiguity of a power expansion in defining the expanded 
function. 

The optimal summation method may be different from whatever we now know. It 
may be Borel, some of its generalizations, Pade, but also something very different. The 
criterion should not be mathematical simplicity or elegance of the method or result, but 
adequateness with which the method in question is able to glean all the physical infor- 
mation that the standard conventional field (and perturbation) theory failed to take into 
account and process. 

Summarizing, one can say that the difficulties of the present quantum chromodynamics 
are, first, of technical nature (technical impossibility of calculating higher-order pertur- 
bative coefficients) and, second, of principal nature (ambiguities of asymptotic power 
expansions, renormalization ambiguities). All of them point to insufficiency of the cur- 
rent state of QCD and have in common the feature that they represent a challenge to 
nonperturbative quantum chromodynamics, without which no essential progress seems to 
be possible. 
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